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ABSTRACT 

Observations indicate the presence of a massive black hole in the Galactic center. Within a half- 
parsec from the Galactic center, there is a population of coeval young stars which appear to reside 
in a coherent disk. Surrounding this dynamically-cool stellar system, there is a population of stars 
with a similar age and much larger eccentricities and inclinations relative to the disk. We propose a 
hypothesis for the origin of this dynamical dichotomy. Without specifying any specific mechanism, 
we consider the possibility that both stellar populations were formed within a disk some 6±2 Myr 
ago. But this orderly structure was dynamically perturbed outside-in by an intruding object with 
a mass ~ W^Mq, which may be an intermediate-mass black hole or a dark stellar cluster hosting 
an intermediate-mass black hole. We suggest that the perturber migrated inward to ~ 0.15 — 0.3 pc 
from the Galactic center as a consequence of orbital decay under the action of dynamical friction. 
Along the way, it captured many stars in the outer disk region into its mean-motion resonance, forced 
them to migrate with it, closely encountered with them, and induced the growth of their eccentricity 
and inclination. But stars in the inner regions of the disk retain their initial coplanar structure. 
Quantitatively, a perturber on a low-inclination or overhead orbit to the disk plane can reproduce 
the observed kinematic structure of these young stars. But this process is unlikely to produce the 
controversial two-disk structure. We predict that some of the inclined and eccentric stars surrounding 
the disk may have similar Galactocentric semimajor axis. Future precision determination of their 
kinematic distribution of these stars will not only provide a test for this hypothesis but also evidences 
for the presence of an intermediate-mass black hole or a dark cluster at the immediate proximity of 
the massive black hole at the Galactic center. 

Subject headings: Black hole physics: Galaxy: center-stars: kinematics and dynamics: stellar dynamics 



1. INTRODUCTION 

Recent observations show that a large number of early- 
type stars, either 0/W-R stars or less massive B stars 
(totally about 90), reside within the central half-parsec 
region of t he Ga l actic center (GC) (iGenzel et al.l l2003l: 
iGhez et all I2003L l2005l : iPaumard et al.ll2006l . and refer- 
ences therein) . These young stars have some distinc- 
tive features (iLevin fc Beloborodovl 120031 : iGenzel et al.l 
l2Q03HPaumard et al.ll2006i r (a) About 26 stars are found 
to reside in a well-defined and moderately thin disk, 
which rotates clockwise in projection and has a scale 
height-to-radius ratio h/r ^ 0.12 ± 0.03. According to 
IPaumard et al.l ()2006f ). about a dozen stars reside in a 
less well-defined counterclockwise rotating disk and these 
two disks are orientated at a large angle (~ 110°) with 
each other. Note that the existence of the counterclock- 
wise rotating disk is still in controversy; nevertheless, ob- 
servations indicate that the motion of some stars (prob- 
ably located at the outer region, see item c below) is 
non-coplanar with respect to that of stars in the clock- 
wise rotating disk, (b) The stars in the clockwise ro- 
tating disk are on low-eccentricity orbits (with eccen- 
tricity e 0.2 — 0.4 on average). In contrast, many 
stars in the counterclockwise disk appear to be on high- 
eccentricity orbits (say, e ~ 0.8). (c) The clockwise ro- 
tating disk is compact and has a sharp inner edge at 
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~ 1" (~ 0.04 pc with an assumed GC's heliocentric dis- 
tance of 8kpc), while the counterclockwise component 
has a ring-like structure which is located further out at 
~ 4" (~ 0.15 pc). (d) The stars in the "two disks" are 
coeval with an age of 6 ± 2 Myr and probably formed 
within a time span of < 2 Myr. (e) A cluster of main- 
sequence B stars (4M0 < < 15 Mq) or the so-called 
"S-cluster" exists in the inner 0.04 pc region, which 
has a spatially isotropic distribution and is distinct from 
the more massive sub-parsec (~ 0.04 — 0.5 pc) young 
stars. Some of these S-stars have high eccentricities (with 
e up to 0.9-1) while others have moderate eccentrici- 
ties. The S-star orbits have provided strong evidence for 
the existence of a massive black ho le (MBH; with mass 
Mbh ^ 3.6 X 106 M,.,) in the GC (iSchodel et al.l l200l : 
iGhez et aLll2005l : lEisenhauer et al.ll2005l . and references 
therein). The young stars in the GC, together with the 
central MBH, provide an interesting dynamical system to 
study. Specifically, the dynamical architecture of the GC 
resembles that of the solar system, which is composed of 
a thin disk of major planets and a thick population of 
minor planets orbiting around the Sun. As studies of the 
dynamics in our solar system have provided us consider- 
able insights on the formation and evolution of the Sun 
and its surrounding planets, investigations on the dy- 
namical system in the GC may also provide us insights 
into the structure, and the formation and evolution of 
the central MBH and stars in the nucleus of our own 
Galaxy and further in general galactic nuclei. 
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The coevality of the sub-parsec young stars suggests 
that they have a common origin. Since the tidal force of 
the central MBH may prevent formation of young stars 
through the collapse of self-gravitating c old molecular 
gas clouds in the vicinity of the MBH (e.g- JSanderslllOO^ : 
lMorrisl[l993( ). other hypothesis have been proposed to 
explain the formation and youth of the stars, including 
some non-conventional in situ formation scenario and 
in-spiraling young star cluster scenario. The in situ 
formation scenario is based on the assumption that the 
young stars were formed via the onset of gravitational 
instability and fragmentation in a hypothetical accretion 
disk around the MBH which no longer exist today. With 
a sufficiently large gas surface density, the self-gravity 
of the perturbation may o vercome the im pedim ent of 
the MBH's t i dal fo r ce fe.g..lLeyin fc Belo borodo^ l2003l : 
iGenzel et all [200l iGoodmanI 120031 : iLevinI 120071) . In 
the in-spiraling star cluster scenario, it is assumed that 
the young stars were originally formed in a dense star 
cluster outside the central half-parsec and transported 
to their present locati on by the effect of dynamical 
friction ()Gerharcil |2001D . An intermediate-mass black 
hole (IMBH; a few thousand solar mass) may be required 
to stabilize the cluster core against its tidal disruption 
by the MBH before they move t o the ce ntral half-parsec 
region (iHansen fc Milosavlievid 2003; Kim et all l2004l: 
Giirkan fc Rasid l2005l : iMcMillan fc Portegies ZwartI 
20031) . Note that the existence of stellar disks in 
observations itself is not sufficient to differentiate these 
competing scenarios since a young stellar disk may form 
in the vicinity of the MBH in both of the scenarios. (For 
a review and more discussions on the p ros and cons of 
these scenarios, see ffmizcl ct a l. 2003; Ghez et al.ll2005l : 
iPaumard et"aLl[2006l ; iAlexandeill2005l .) 

Besides their young age, the peculiar orbital distribu- 
tion of these stars also provides clues to their formation 
mechanism and the dynamical structure in the GC. It 
can also be used to differentiate various models. For 
example, in the in-spiraling star cluster scenario, the in- 
teractions between the central I MBH and the stri pped 
young stars have be e n stu died bv lLevin et al.l (|2005| ) and 
iBerukoff fc HansenI ()2006D . Although they were able to 
simulate the orbital features of stars on the thin disk, 
they had difficulties to account for the kinematic prop- 
erties of those stars with eccentric orbits in the coun- 
terclockwise rotating disk. In this paper, we propose a 
dynamical model to account for the orbital distribution 
of the sub-parsec young stars in the GC. Rather than si- 
multaneously investigating the formation mechanism of 
these young stars and their orbital distribution, we as- 
sume that all the sub-parsec young stars were formed in 
a clockwise rotating disk initially (e.g., due to instability 
or fragmentation of a massive accretion disk) with small 
eccentricities. We show in § [3] that, an isolated stellar 
disk by itself is unlikely to evolve from the presumed 
circular orbits in a thin disk to the currently observed 
multi-component orbital configuration. This inference 
is based on the determination that the timescale (e.g., 
the two-body relaxation timescale and the resonant re- 
laxation timescale) for a single stellar disk to relax with 
respect to the background stars is usually substantially 
longer than the age of the young stars. Thus, the life-long 
eccentricity growth of the disk stars due to interactions 
among the stars is insignificant (not substantially larger 



than 0.2; [Alexander et all [20071 ). Therefore, the exis- 
tence of young stars with high eccentricities and high in- 
clination angles relative to the inner stellar disk requires 
an explanation, at least in the disk formation scenario of 
the young stars in the GC. 

It is entirely plausible that, under the action of dynam- 
ical friction, a (dark) star cluster and/or an IMBH oc- 
casionally sink into the proximity of the GC. The (dark) 
star cluster (which may also have a central IMBH) or 
the IMBH plays the role of a massive perturber (or in- 
truder), spiraling inward, inducing asymmetric perturba- 
tion of the gravitational potential and exerting torques 
on the disk stars. Interactions between the perturber and 
the disk stars, through either close encounters or some 
resonant interactions, may change the orbital distribu- 
tion of the young stars. (In contrast to the in-spiraling 
young cluster hypothesis, the primary role of the cluster 
is to dynamically excite the eccentricity and inclination 
of the stars along its path rather than to directly deliver 
its own young stars to the GC.) In this paper, we develop 
dynamical models in order to examine whether the single 
stellar disk initially formed may evolve into the current 
observed orbital configuration, especially the eccentricity 
and inclination distribution, during and after the passage 
of an inward migrating massive perturber. 

This paper is organized as follows. In § [21 we review 
the components observed in the GC. In §[3l we list some 
relevant dynamical timescales in the GC. In § [31 we de- 
scribe dynamical models of gravitational interaction be- 
tween the young stars in a primary disk and a massive 
perturber migrating inward from outside the stellar disk 
in the GC. We present the results of a set of numerical 
simulations of the dynamical models. Then we compare 
the model results with observations. Finally, conclusions 
are given in § [5l 

2. COMPONENTS IN THE GC 

We consider three components which contribute to the 
gravitational potential near the center of the Galaxy; a 
central MBH, a stellar cusp of old stars, and a population 
of young stars mentioned in § [l] 

The MBH has a mass A/bh ^ 3.6xl0'^Mo (|Ghez et al.l 
[2005t [Eisenhauer et al.[|2005D . The radius of the central 
MBH's sphere of influence is defined by 

GA/bh , ^ /100kms-i\' 
an = — — ~1.6pcmBHl I , (1) 

where G is the gravitational constant, a is the one- 
dimensional velocity dispersion, and ttibh ~ A/bh/(3.6 x 
IO^Mq) in the GC. 

The stellar cusp surrounding the young stars is mainly 
composed of old stars and stellar remnants including 
stellar-mass black holes and neutron stars. The mass 
density of th e stellar cusp in the GC can be described by 
a power law (jCenzel et al.ll2003h . i.e., 

p(r) = po(^^) , 1^1,2, (2) 

where po = 1.2xlO'^M0 pc'^, = 0.4pc, ai = 1.4±0.1 
for r < ro and a2 = 2.0 ±0.1 for r > vq. This mass 
density gives the enclosed stellar mass in the inner cusp 

by 

M.«.) = l^(^)"" 
3-ai \roJ 
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and 



MJ< r) 



= eA/BH/i(a;), if r < tq 

3 - ai 3 - a2 
:eMBH/2(a;), if r > tq 



(3) 



(4) 



where e = 



and f2{x) 



^/A/bh ^ 0.16, X - r/ro, /i(x) 



^,3 — ai 



3 — ai 



. The gravitational po- 



3 — £12 

tential due to the cusp stars is 

^ , , GAM<r) , (3 - ai)eGMBH , 

$*(rj = 1 in 

r To Vro 



for r > ro and 

GAf,(<r) (3-ai)eGA/BH 



(2-ai)ro 



, (5) 



(6) 

for r < tq. The total potential due to the MBH and the 
cusp stars is $(r) = — + <i>*(r). Since the young 

stars interested in this paper mainly reside within the 
central half parsec and tq = 0.4 pc, for simplicity, we 
adopt a2 = «! ~ 1.4 in our calculations below and use 
equations Q and ([5|) even for r > tq. Our main results 
will not be affected by slightly changing the value of the 
cusp slope. 

3. SOME RELEVANT DYNAMICAL TIMESCALES IN THE 

GC 

In this section, we estimate some dynamical timescales 
in the GC, which may be used for references or to justify 
the approximations in the dynamical models described 
in §11 

• The orbital timescale of a star rotating around the 
central MBH is 



GMbh 



1 r in3 -1/2 

1.5 X 10 yr 



1.5 



(7) 



^O.lpc 

where a is the orbital semimajor axis of the star. 

• The local two-body relaxation timescale Treiax is 
given by 



0.34cr3 



' relax ' 



G^p* (m*) In A 
2.0 X lO^yr 



In A 



^BH ' 



BH 



IM,. 



3/2 



M4< r) (m*) V0.4pc 

q 1/9 10 IMpj ' - ^ ""-^ 
'1.2 X lO^yr mB'j:f " 



-In A 



0.4pc 



(8) 



where (m*) is the me an stellar mass and p,, is 
the stellar densi ty (e.g.. iBinnev fc TremaindriOSTl : 
lAlexanderlll999| ). In the estimate of this timescale 
below (in Fig. [T]), we set (m*) ~ 1 Mq and InA ~ 
10. 



• The precession due to the stellar cusp: the stellar 
cusp (mentioned in § [2]) may cause stars deviated 
from purely Keplerian motions and introduce apsi- 
dal precession of their orbits. This apsidal preces- 
sion timescale can be estimated by 



TnCUSp 

prcc ' 



Af, 



BH 



Af,(< r) 



T'orb(r-) 



'7.5 X lO^yr va^^ 



0.4pc 

with assuming Af,(< r) A/bh- 



-0.1 



(9) 



• Resonant relaxatio n timescal es (see details in 
Rauch fc Tremain^ [l99 6: H opman fc Alexander! 



20061) : in a timescale much longer than the or- 



bital period and shorter than the apsidal precession 
timescale mentioned above, the star can be approx- 
imated by a fixed wire whose mass is the stellar 
mass, whose shape is a Keplerian ellipse, and whose 
linear density is inversely proportional to the local 
speed in the elliptical orbit. The wires precess on 
the timescale TpJJfP and exert mutual torques which 
induce angular momentum relaxation. The cumu- 
lative effects of the torques result in the change of 
the absolute value of the angular momentum by a 
timescale of 



- 3.56 X 



A/, 



BH 



M4< r) 



(to*) 



Torb(r) (10) 



(hereafter the scalar-resonant relaxation timescale) 
or the change of the direction of the angular mo- 
mentum vector by a timescale of 



^rcs.v . ^ A/BH + Af,(<r) ro,b(r) 

^rclax -^U.O^X ^^^^^ iVl/2(<r) 



(11) 



(hereafter the vector-resonant relaxation 
timescale), where A^(< r) = M^,{< r)/ {m^). 

The sub-parsec massive young stars is a different 
population from the major component (old stars 
with (m,) ~ IM0 set here) in the Galactic center. 
We still take the above equations as the resonant 
relaxation timescales of the young stars, since their 
number and total mass are much smaller than the 
background old stars and the effect of the resonant 
relaxation on the young stars should mainly be due 
to the old stellar populat ion (see the derivation in 
iRauch fc Tremainel[l996l ). 

• The apsidal precession due to the general relativity 
correctio n to the Newtonian equations of motion is 
given by (jMisner et al.lll973D : 

ac^{\ — e^) 

3GA/BH 

= 3.0xl0VrTO3^/^(l-e^) (^)'''xi2) 

where c is the speed of light. 

• The Lcnsc-Thirring precession timescale is given by 
iMisner et al. 1973|) 



Tut — 



2^3(1 _ e2)3/2 



2J 
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.2,2xlO'»y,»-Ve5(l-.V'^(i^) . 

(13) 

where J is the angular momentum of the MBH 
and a* = J/{GM^-^/c) is the dimensionless spin 
parameter of the MBH. 

The above timescalcs are plotted as a function of the 
distance from the central MBH in Figure [TJ with us- 
ing the parameters of the observe d stellar cusp described 
in § [2] (see a similar diagram in iHopman fc Alexander! 
|2006[ ). As shown in Figure[TJ the timescales Trdax, ^reiax ' 
and Tjfgjax' ^'"^ much longer than the stellar age Tage 
at r > 0.04 pc, which suggests that the spatial distri- 
bution of the original disk located at 0.04 — 0.5 pc can- 
not be erased through either two-body or resonant relax- 
ation processes. The apsidal precession timescale Tgr 
is much longer than Tagc for the stars at a > 0.1 pc 
and it is comparable to the age of the disk stars with 
moderate eccentricity at a < 0.04 pc. The apsidal 
precession may affect the orbits of the S-stars (since 
Tqr < T'ago for the S-stars), but it should be insignifi- 
cant to the orbital evolution of the disk stars located at 
a ~ 0.04 — 0.5 pc. The Lense-Thirring precession is gen- 
erally not important for those disk stars since Tlt ^ Tagc 
at a 0.04 — 0.5 pc, b ut it may be important for th e in- 
nermost two stars S2 (jLevin fc Beloborodwll2003f l and 
S14 (Tlt - 4.1 x lO^aji yr for S2 ~ 8.9 x lO^aJ^ yr for 
S14). As seen from this Figure, the apsidal precession 
timescale due to the stellar cusp is much shorter than 
the age of the stars but much longer than their orbital 
period. This precession affects the secular orbital evolu- 
tion of the young stars in the stellar disk. 

According to the estimates of these relevant timescales 
(see Fig. [TJ, the disk stars in the GC may be described 
to be moving in a stellar cusp with a smooth potential 
(the apsidal precession due to the stellar cusp is nat- 
urally included in the description), and the effects of 
other precessions, apsidal precession due to the general 
relativity correction and the Lense-Thirring precession, 
may be neglected in the following calculations. In this 
paper, we do not pursue a study on the origin of the 
kinematics of the S-stars (with age of a few 10^ yr), for 
which the resonant relaxation, as well as the precession 
due to the general relativity correction and the Lense- 
Thirring preces sion, arc invo lved (see als o discussions in 
iHopman fc Alex ander 2006i; iLevinI l2007f ) and a smooth 
potential may not be a good approximation any more. 

In addition, if the young stars were formed in a gas 
disk, we assume that the gas disk was depleted quickly 
and gas drag is not important for their orbital evolution 
(damping the orbital eccentricity) at least for the two 
reasons based on observations: (1) most of the stars were 
formed in a short timescale < 2 Myr compared to their 
age 6±2 Myr; and (2) the MBH is accreting material via 
a rate around 1 0~^ — 10~^ Mm yr~^ through a tenuous 
thick disk fe.g.. lMelia fc Falckell200lt ). which should not 
have any significant drag on the motion of the young 
stars. 

4. MODEL AND NUMERICAL EXPERIMENTS 

It is plausible that occasionally there may be a star 
cluster (and/or an IMBH) in-spiraling into the central 
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Fig. 1. — Relevant timescales for the young stars in the GC as 
a function of Galactocentric raidus r or semimajor axis of a stellar 
orbit a: the two-body relaxation timescale T^eiax (blue solid line), 
the scalar-resonant relaxation timescale T^^^^^ (red solid line), the 

vector-resonant relaxation timescale T^^i^ (magenta solid line), 
the Lense-Thirring precession timescale Tlx by assuming the di- 
mensionless spin parameter of the MBH a« = 0.5 and the eccen- 
tricity of the stars e = 0.5 (cyan solid line), the apsidal precession 
timescale due to the general relativity correction Tqji with as- 
suming e = 0.5 (green solid line); the apsidal procession timescale 
due to the stellar cusp Tp^J? (dashed line); and the orbital pe- 
riods of the young stars T^rh (dotted line). For comparison, the 
age of those young stars Tagc ~ 6 it 2 Myr and the inner radius 
of the clockwise disk rinncr ~ l" ~ 0.04 pc arc marked by the 
dot-das hed lines. For reference, the rele vant timescales for several 
S-stars l lEisenhauer et all 120051 : IGhez ct al., ,20051) with measured 
orbital parameters are also marked in this Figure (filled black cir- 
cles: orbital periods; cyan triangles: Lense-Thirring timescales; 
green squares: apsidal precession timescales due to the GR correc- 
tion). See detailed discussions in §|3] 

regio n of the Galaxy (see also discussions in lPerets et al.l 
l2007ft . In this section, we study how the orbital config- 
uration of the young stellar disk are affected by such a 
perturber. To isolate the problem, we first estimate how 
the orbits of the young stars is affected by a perturber 
rotating around the MBH at a fixed distance in §STJ We 
consider a range of inclination (from 0° to 180°) between 
the perturber's orbital plane and the stellar disk plane. 
In § 14. 2[ we illustrate how the orbital configuration of 
the stellar disk is affected by the inward-migration of a 
perturber. 

4.1. Perturbation on the stellar orbits due to a massive 
perturber at a fixed distance from the central MBH 

The orbital precession caused by the torque exerted by 
a massive perturber at a fixed distance from the central 
MBH (i.e., on a circular orbit) changes the inclination 
angles of the young stars, but not their eccentricities. 
In order for the clockwise rotating young stars in the 
inner region to retain their nascent disk's initial configu- 
ration, the oscillation amplitude of the inclination angle, 
induced by the nodal precession, needs to be small (say, 
< 20°; this number is roughly the maximum range of the 
estimated inclination angles of the yo ung stars in the in- 
ner d isk relative to the disk normal; iBeloborodov et all 
l200(]( l. The oscillation amplitude should be small if the 
inclination angle between the orbital plane of the mas- 
sive perturber and the stellar disk, denoted by /3, is close 
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to 0°or 180°, but it can be large for other intermediate 
values of /3 unless the precession timcscalc is much longer 
than the stellar age. 

Approximating the gravity of the perturber and the 
stars as ring s, their nodal precession frequency is e sti- 
matcd to be ()Navakshinll2005t iNavakshin et al.(l2OO60 



(r2 + r2)3/2 r 



(14) 



where Mp is the mass of the perturber which is much 
more massive than a young star, Qk is the Kep- 
lerian angular frequency of the stars, S ~ 



sm sm ( 



27r [1 — (5(cos (3 sin tp' sin (/)+cos 0' cos 0)]'^/^ ' 

and I{S,f3) = I{S,tt — (3). Using equation (O, we 



show the precession timescale Tprocoss = 1/w* oc 1/Mp 
in Figure [2] with different parameters and /3. As seen 
from Figure [21 for a perturber with mass Mp = 10^ Mq 
and located at a relatively large Galactocentric distance 
(rp ~ 0.6 pc; blue Hues), the precession timescales in 
some region (r ~ 0.5 pc) may be smaller than the stellar 
age for some values of (3 and the orbital inclinations of 
the stars in this region may have substantial variations 
due to the nodal precession; but the precession timescale 
at r ^ 0.04 — 0.15 pc is generally much longer than the 
stellar age and the initial coplanar structure of the inner 
disk can be well preserved for any values of /3. If the 
perturber is located at a smaller Galactocentric distance 
(rp ~ 0.2 pc; red lines), the nodal precession timescale 
of the young stars at r ^ 0.1 pc is much shorter than 
or comparable with their age if /? is in the range ~10°- 
80° (or ~100°-170°); thus the orbital incHnations of the 
stars not only in the outer region ~0.15-0.5 pc but also in 
the inner disk may have substantial variations or a warp- 
ness may be developed in the stellar disk due to the nodal 
precession (these qualitative estimates are confirmed by 
direct numerical calculations of the stellar orbital motion 
in ? l4.2p . If the perturber's orbit is almost perpendicular 
to the initial disk plane (red long-dashed line), the pre- 
cession timcscalc for the young stars in the inner region 
is long enough compared to their age, and thus the inner 
disk may be not affected by the precession and its ini- 
tial coplanar structure can be preserved. If the perturber 
orbit is nearly parallel to the disk plane, as mentioned 
in the paragraph above, although Figure [2] shows that 
the precession timescale is short compared to the stellar 
age, the disk configuration in the inner region may still 
be preserved since the oscillation of the normal to the 
stellar orbital planes around the normal to the perturber 
orbital plane is small. Thus, the preservation of the in- 
ner disk suggests that if the perturber with mass 10"* Mq 
is located at rp ~ 0.2 pc, its orbit should be more likely 
to be parallel or perpendicular to the assumed primary 
disk. This result is also demonstrated in § 14.21 below. 
For a substantially smaller A/p, the precession timescale 
Tprecess Can bc long enough so that the inner disk may be 
maintained for a large range of (3. 

4.2. Perturbation on the stellar orbits by an 
inward-migrating perturber 

In this subsection, we study how the background stellar 
orbits evolve as a massive nearby perturber migrating 
inward. For simplicity, we assume the perturber to be 
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Fig. 2. — The timescales of the nodal precession of stellar orbits 
due to the torque exerted by a massive perturber with mass Mp = 
10* Mq at a fixed distance from the central MBH rp. The stars 
are on nearly circular orbits with the Galactocentric distance r. 
The left (red) and right (blue) sets of lines are for a perturber 
located at rp = 0.2 pc and rp = 0.6 pc, respectively. For each set, 
diff'erent line typos represent different inclination angles /3 between 
the orbital plane of the massive perturber and that of the stellar 
disk [from top to bottom: /3 = 85° (long-dashed linos), 75° (short- 
dashed lines), 45° (dotted lines), and 15° (solid lines)]. The age and 
position of the inner stellar disk (rage ~ 6 it 2Myr and r ~ 0.04— 
0.15pc) observed in the GC are marked as solid green lines. The 
dotted green line illustrates the extension of the assumed original 
disk to ~ 0.5 pc. To preserve the coplanar structure of the stars 
in the region 0.04-0.15 pc not to be significantly affected by the 
nodal precession, the perturber should at least (i) be located at a 



far distance (e.g., 



' 0.6 pc), or (ii) have an orbit nearly parallel 



or perpendicular to the disk plane if the perturber is located at 
a small distance (e.g., rp ~ 0.2 pc), or (iii) have a substantially 
smaller mass. The orbits of the stars at the outer region (r ~0.15— 
0.5 pc) can generally be affected by the nodal precession at some 
values of /3 unless the perturber mass is substantially smaller. See 
detailed discussion in ij 14.11 

an IMBH with a point mass potential. If the perturber 
is a dark cluster, the gradually intensifying Galactic tidal 
effect is likely to induce the perturber to lose mass during 
its inward migration. We also consider the possibility 
that the perturber has a declining mass. 

Before the perturber becomes bound to the central 
MBH, the perturber spirals into the Galactic center un- 
der the action of dynamical friction which induces a de- 
celeration 



~dt 



V$(r), 



(15) 



where the dynamical-friction timescale 



T = tdf = 



87rG2 lnAMpp(r)[erf(X) - ^exp(-X2)] 



(16) 

^^p = l^^pl is the velocity of the perturber, X = 
^^^E— , and In A is the logarithm of the ratio of 

the maximum an d minimum impact para meters and 
A ~ Mbh/Mp (|Binnev fc Tremaind Il987t) . As the 
perturber or the IMBH migrates inward and forms 
a bound binary black hole (BBH) with the central 
MBH at Tp ~ flH (see eq. [1]), it continues to lose en- 
ergy and angular momentum through dynamical fric- 
tion. But, the influence of dynamical friction on the 
IMBH's orbit becomes less efficient as its orbital pe- 
riod decreases and its orbital velocity increases. Af- 
ter the BBH becomes hard at ah = (A'/p/4AjfBH)aH — 
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0.004 PC mBH(^^^^^) , it loses energy 

mainly through three-body interactions with stars pass- 
ing by its vicinity. The orbital decay timescale of a hard 
BBH in t he GC is about th = rp/fp ~ 6 x 10^ yr (see 
eq. 38 in lYu fc Tremaind l2003f ) . The gravitational ra- 
diation timescale of the BBH is longer th an the Hubble 
time (see eq. 39 in lYu &: Tremaind |2003[) and it is un- 
likely to be significant in the spatial range of the BBH 
considered in this paper. During the transition stage 
(after the BBH becomes bound but before it becomes 
hard), the BBH's hardening timescale may be higher 
than the estimate from the dynamical friction timescale 
tdt(~ 10^ yr at 1 pc for Afp ^ 10^ Mq), as this process be- 
comes less efficient, and it increases to that at the hard 
stage as the IMBH migrates in fYj||2002[) . Under 
the constraint placed on the distance between Sgr A* 
a nd the center of mass of the BBH (eq. 35 or Fig. 2 
in lYu fc Tremain"3l2003[ ). the mass of the secondary BH 
(or star cluster) Mp should be smaller than ^ 0.03Mbh 
if it is located at 0.3 pc and smaller than ^ 0.08A/bh 
if it is located at O.lpc. If the GC (or the central 
MBH) was formed through the assemblage of sinking star 
clusters (or IMBHs, which are considered to be candi- 
dates for the massive perturber in this paper), we can 
also infer their masses on the assumption that we are 
not living in a unique time. Such an argument would 
lead to a mass estimate for a typical perturber to be 
~ [Mbh + M*(< Ipc)] X 10Vr/10^°yr ~ 5 x IO^Mq. 
Finally, star clusters with mass around 10'' Mq are found 
within several tens pc from the GC. Based on these con- 
siderations, we adopt Mp = 10^ Mq as a fiducial mass 
for the massive perturber in general. We also consider 
the potential implications for other values of Mp. 

In the first set of simulations of stellar responses to 
a migrating perturber, its initial orbit is assumed to 
be nearly circular with a radius 1 pc. The path of in- 
ward migration follows equation (|16p. where the mi- 
gration timescale r is set to several different values 
(10^, 10^, 10^ yr). The initial conditions of the perturber 
are set so that it attains an (almost constant) eccen- 
tricity ^ 0.02 after it moves sufficiently close to the 
MBH (rp < 0.2 pc). (In subsequent models, the poten- 
tial implications for the higher values of the perturber's 
eccentricity are also considered. But we do not con- 
sider a highly eccentric orbit for the perturber, since 
the orbital decay of a dark cluster or a secondary BH 
via dynamical friction and three-body interactions with 
the central MBH and stars generally does not result in 
a highly eccentric orbit, e.g., with e ccentricity < 0.3; 
iPohiarev fc Reeslll994l : lQuinlanlll996i ) The orbital evo- 
lution of the perturber is shown in Figure [3l These re- 
sults indicate that 15Myr (or lOMyr) would be required 
for the perturber to migrate inward from 1 pc (or 0.5 pc) 
to O.lpc if T = 10^ yr. In this paper, we do not ad- 
dress the interesting issue whether the migration of the 
perturber may be correlated with the formation of the 
young stars. We simply assume that the young stars 
formed throughout the disk prior to or during the mi- 
gration of the perturber. Note that the time i = does 
not necessarily represent the time when the young stars 
formed. 

We simulate the orbital evolution of 75 test particles 
as young stars, with the logarithm of their initial semi- 




10 20 30 40 50 

I (Myr) 

Fig. 3. — The orbital evolution of a perturber used in the dy- 
namical model proposed in § 14.21 The perturber is initially on a 
nearly circular orbit at 1 pc and its evolution follows equation l|15| l 
with T = 10*, 10^, 10*^ yr from top to bottom. 

major axes a distributed uniformly in the radial range 
lg(0.04 pc)-lg(0.5 pc). Each star moves independently in 
the potential of the central MBH, the inward-migrating 
IMBH, and the stellar cusp. We set the inclination of the 
perturber to zero and it does not change in the spherical 
gravitational potential of the cusp. Five sets of initial 
inclination angles of young stars i are chosen randomly 
in the range 0°-10°, 20°-30°, 80°-90°, 110°-120°, and 
170°-180°, respectively. Their other orbital elements rel- 
ative to the central MBH, such as the longitude of as- 
cending node, argument of pericenter, and true anomaly, 
are chosen randomly in the range [0° , 360°] . The initial 
eccentricities e's of the orbits around the MBH are set 
to zero. But, they rapidly attain finite amplitudes, due 
to the gravity contributed by the stellar cusp. Given the 
position and velocity of a test particle at a given time, we 
can still define its time (or position)-depcndent semima- 
jor axis and eccentricity by assuming that the particle is 
on an elliptical orbit around the MBH at the appropriate 
time. In order to illustrate the change of the orbital con- 
figuration of a test particle, the angle of the orbital plane 
of each particle relative to its initial plane, 0, is used be- 
low, which may describe the observed inclination angles 
of the young stars relative to the observed inner disk (our 
conclusions do not change if the initial longitudes of as- 
cending node of the particles in the simulations are set to 
be the same so that all the particles are initially on a thin 
disk). During their evolution, the orbits of two test par- 
ticles with the same inclination angles do not generally 
lie on the same plane. 

We first show the simulation results for a model in 
which the perturber lies almost on the initial orbital 
plane of the disk stars (z S [0°,10°]). We set the per- 
turber mass to be A/p = 10^ M© and the migration 
timescale to be r = 10^ yr. From the simulation re- 
sults, we find a range of stellar responses. For illustration 
purpose, it is useful to classify these responses into the 
following categories with representative cases shown in 
Figures [4H3 

A) Figure [4] shows the evolution of a test particle 
(i S [0°, 10°]) which is captured into the 2:1 mean-motion 
resonance of inwardly migrating perturber, i.e., their or- 
bital period ratio becomes 2:1 at t/r ~ 1. After being 
captured, the particle migrates inward with the decay- 
ing perturber. During the resonant march, the particle 
retains an adiabatic invariant which for a (p-l-9,p) mean- 
motion resonance is C = a}^'^[{p + q) — p{l — e^Y^"^ cosz] 
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Fig. 4. — The orbital e volu tion of a test particle in the dynam- 
ical model proposed in § 14.21 The initial inclination angle of the 
particle is close to 0. From top to bottom, the panels show orbital 
period ratio to the perturber, semimajor axis, eccentricity, incli- 
nation angl e to the orbital plane of the perturber, the quantity C 
defined in S I4.2I argument of pericenter, and the angle between the 
orbital plane of the particle and its initial orbital plane 8. The 
particle is first captured into the 2:1 mean-motion resonance and 
then released from the resonance but into the uj = ±90° secular 
resonance. 
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Fig. 5. — The orbital evolution of a test particle that is captured 
into the 2:1 mean-motion resonance and then into ui — Op = 0° 
secular resonance, where ui — Qp is the longitude difference of peri- 
center between the particle and the perturber. 



Fig. 6. — The orbital evolution of a test particle that initially 
follows the evolutionary pattern as that in Figure |4] and then is on 
an irregular orbit. 



(|Yu fc Tremaindl2001h . The magnitude of p > and 
q > —p are integers, and p = q = 1 for the inner (2,1) 
resonance (by "inner" refers to the resonances he inside 
the perturber's orbit, i.e. with a < Op). For this adia- 
batic invariant, a reduction in the particle's semimajor 
axis during the migration can be compensated by an in- 
crease in its e and i. The results in Figure |4] indeed 
indicate a growth in the particle's eccentricity. When 
its eccentricity is sufficiently high (^ 0.9 at t/r ~ 1.5), 
its inclination also grows. A comparison between the 
results of the two sets of simulations indicates that the 
inclusion of the stellar cusp in the gravitational potential 
promotes the elevation of the particle's inclination to a 
range between 0°-90°, rather than to a polar (i = 90°) or 
a retro grade orbits (i = 180° ) around a point-mass po- 
tential |^T5I^emaini[2Q0l|). The asymptotic inclina- 
tion decreases with increasing fractional contribution of 
the stellar cusp mass to the gravitational potential. Dur- 
ing the lifting of its inclination, the particle enters into a 
secular resonance with perturber's argument of the peri- 
center uj, librating around ±90°. At an advanced stage 
of evolution (t/r ^ 1.7), the particle is released from the 
perturber's mean-motion resonance when it attains rela- 
tively large eccentricity and inclination. But, the particle 
retains its secular resonance {to ^ ±90°) with the per- 
turber which prevents close encounters between them. 
Consequently, the orbit of the particle is stable on the 
timcscalc of our integration (^ 5 x 10^ initial orbital pe- 
riods of the perturber) . 

B) Figure [5] illustrates the evolution of second repre- 
sentative particle with a similar i e [0°, 10°] but have a 
different initial location to that in the previous example. 
The initial evolution of this particle is generally the same 
as that in the previous model (Figure (H). The particle is 
captured into the perturber's mean-motion resonance of 
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Fig. 7. — Orbital distribution of simulated test particles when the perturber with mass 10* migrates to 0.15 pc (denoted by the 
subscript "f " ; the initial values of their orbital parameters are denoted by the subscript "i" ) . The initial inclination angles of the particles 
to the perturber plane are in the range [0°, 10°]. Note that the distance to the central MBH »"/ is generally different from the somimajor 
axis of the particle a j due to the non-zero eccentricity e j . Sec § 14.21 



the perturber at t/T>^ 1.1 (see Figured]). At t/r ^ 1.4, 
the test particle enters into a state of secular resonance 
with the perturber such that the difference of their peri- 
center longitude rj = Cj ~ ujp (where <!» = il + and 
is the longitude of ascending node) librating around 0. 
During this phase, the quantity C is no longer an in- 
variant. The capture of the particle into this perturber's 
secular resonance can be understood through the evolu- 
tion of the difference between the precession rates of the 
perturber and the test part icle, which may be described 
by d-q/dt = Ai + A2 003(77) (|Nagasawa. Lin fc Ida, 20031 ). 
The right-hand side of this evolution equation includes 
the contribution from the precession of the test particle 



induced by the secular interaction with the perturber and 
also the contribution from the precession of both the per- 
turber and the test particle induced by the stellar cusp. 
As the perturber and the test particle migrate inward, 
both of these contributions evolve with time. At some 
point during the course of migration, Ai <^ A2 and the 
longitude difference of pericenter librates around 77 = 
(or 180°). The eccentricity ratio of the test particle and 
the perturber e/ep [which evolves as d{e/ep)/dt ex. sin 77; 
iNagasawa. Lin fc Idall20d3j also librates around an equi- 
librium value. As the test particle continues to migrate 
inward with the perturber, it may retain the state of 
secular resonance but the equilibrium value of the ec- 
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Mn = 3 X 10'* Mq migrates to r 

centricity ratio changes as the Hamiltonian of the sys- 
tem evolves. Similar s ecular resonance was discussed in 
iNagasawa. Lin fc Idal (|2003f ). where the depletion of a 
protostellar disk around a planetary system causes the 
change of the precession rates of the system. In this pa- 
per, the change of the precession rates is caused by the 
migration of the system in a stellar cusp. 

C) Figure [6] shows a third example of the evolution of 
a test particle {i e [0°, 10°]), which has a similar evolu- 
tionary pattern as that in Figure [Hat t/r < 1.4 and is 
on an irregular orbit afterwards. 

The above three categories essentially convert the en- 
tire domain of stellar responses. But the actual outcome 
of the perturbcr's passage sensitively depends on the par- 



ticles' initial orbital parameters. In this set of simula- 
tions, most of the particles with initial semimajor axis in 
the range 0.2-0.5 pc fall into case C, and the stars with 
smaller semimajor axes generally fall into case A and a 
few in case B. Figure [7] shows the distribution of the or- 
bital parameters of the test particles. In this figure, the 
subscripts "i" and "f " of various variables represent their 
initial and final values, respectively. For the "final" val- 
ues, we refer to those when the perturber migrates to 
Tp = 0.15 pc {t/r ~ 1.1). (Note that the value of rp 
is an instantaneous position of the perturber, and the 
perturber does not necessarily stop there but continues 
to migrate inward under the action of dynamical fric- 
tion with time going on.) As seen from Figure[7jg), (d). 
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and (h), the scmimajor axes, eccentricities, and inclina- 
tion angles of particles with small initial scmimajor axis 
(oi < 0.1 pc) arc not greatly affected by the migration 
of the perturber. Test particles with initial scmimajor 
axes in the range 0. 1-0.3 pc have essentially the same fi- 
nal scmimajor axes (~ 0.1 pc) because these particles are 
captured into the 2:1 mean-motion resonance and forced 
to migrate with the perturber during its orbital evolu- 
tion. As seen from panels (d)-(g), their eccentricities 
and inclinations may be excited to high values. Panel 
(c) indicates that the test particles with high inclination 
angles generally have high eccentricities. The particles 
with the same scmimajor axes may have different instan- 
taneous distances r/ from the central MBH due to the 
non-zero eccentricity e/. Panels (a) and (b) show the 
distribution of the eccentricities and inclination angles 
versus the distances. Panel (a) indicates the tendency 
that the eccentricities of the test particles are low at the 
small radii and high at the large radii. Panel (b) indicates 
that the stellar disk may be maintained at the inner radii 
[lg(r//pc) < —0.8] and the test particles at the outer 
radii have high inclination angles relative to the inner 
disk. Most particles with initial scmimajor axes larger 
than ~ 0.35 pc are captured into 2:1 resonance during 
the early stages of the perturber's evolution but they are 
released from the resonance and attained irregular orbits 
after the perturber has migrated inside rp = 0.15 pc. 

We also check the dependence of the final orbital dis- 
tribution on the perturber's position rp, mass Mp, and 
eccentricity. For a smaller rp = 0.1 pc, the inclination 
angles and eccentricities of the stars at the inner radii 
(r/ < 0.1 pc) may also be excited to high values and the 
inner disk can no longer be maintained. For a larger 
rp ~ 0.3 pc, although the 'final' eccentricity distribution 
is qualitatively similar to those shown in Figure [71 but 
the inclination angles 9f cannot be excited beyond 20°. 
A relatively low-mass perturber (e.g., Mp = 10^ Mq) 
cannot significantly excite the inclination angles, either, 
although it can excite some particles' eccentricities to 
high values by capturing them into the 2:1 or 3:2 mean- 
motion resonance. An increase in the perturber mass 
(e.g., Mp = 3 X 10^ Mq) enhances its secular perturba- 
tion and widens its mean-motion resonances. So does 
an increase in the perturber's eccentricity or an increase 
in the test particles' initial eccentricities of the test par- 
ticles. In these cases, some particles are captured into 
mean-motion resonances other than the 2:1 resonance, 
such as, the 3:2, 3:1, or 4:1 resonance. Some particles 
may transit from one mean-motion resonance to another 
mean-motion resonance during their orbital evolution. 
Relatively more particles are captured into secular reso- 
nance (D — (Dp = 0° or 180°. 

We now consider the possibility that the perturber mi- 
grates to the proximity of the stellar disk on a highly in- 
clined orbit. We show the case that the initial disk plane 
is almost perpendicular to the perturber orbital plane in 
FigurelH for which Mp = 3 x 10'' Mq and rp = 0.3 pc. In 
this case, most of the main features of the distribution 
in Figure [7] are preserved, including the maintenance of 
a coherent inner stellar disk with Of < 20°. The angles 
of the final orbital plane of the stars relative to their ini- 
tial disk plane 6f have a much larger range 0°-180° than 
those shown in Figure [71 Figures M and [TUl illustrate 
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Fig. 9. — The orbital evolution of a test particle with initial 
inclination angle close to 90°. The particle is first captured into 
the 2:1 mean-motion resonance and then is on an irregular orbit. 
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Fig. 10. — The orbital evolution of a test particle with initial 
inclination angle close to 90°. Although the inclination angle of 
the particle relative to the perturber orbital plane i changes little, 
the angle of the particle orbital plane to its initial orbital plane 
6 changes significantly due to the nodal precession in the gravita- 
tional potential of the stellar cusp and the change of the inclination 
angle from an initial value close to 90° to ~ 70° at t ~ 0.75t caused 
by a close encounter. The angle 8 is almost constant or does not 
oscillate after t/r ~ 1.5 because the perturber has migrated into 
the inner region and does not have significant effects on the orbital 
evolution of the outside particle anymore. 
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Fig. 11. — Orbital distribution of simulated test particles with initial inclination angles i G [20°, 30°] when the perturber with mass 
10* Mq migrates to 0.15 pc. Other initial conditions are the same as those in Figure [7] 



two representative particles' orbital evolution. In Fig- 
ure [9l after the particle is captured into the perturber's 
2:1 mean-motion resonance (or 3:2 resonance for some 
other particles), both its eccentricity and its inclination 
also increase as the particle migrates inward. The mag- 
nitude of C is an adiabatic invariant at t/r < 0.8 and 
the particle is on an irregular orbit after it undergoes 
a close encounter with the perturber at that epoch. In 
Figure [TOl the inclination angle of this particle i changes 
little during the evolution; but after i flips from an ini- 
tial value close to 90° to ~ 70° because of its close en- 
counter with the perturber at t/r ~0.7-0.8, the angle 
between the orbital plane of the particle and its initial 
orbital plane changes significantly due to the nodal pre- 



cession of the particle orbital plane. The excitation of 
the eccentricities of the test particles shown in Figure [8] 
results from their capture into mean-motion resonances 
and close encounters with the perturber. In addition to 
resonance capture and close encounters, nodal precession 
due to the perturber also results in the excitation of the 
inclination angles 9. The distribution of the final semi- 
major axes of the particles in the outer region (e.g., see 
Qf > 0.15 pc in Fig. [8^) is broader than that shown in 
Figure [7] due to the scattering through close encounters 
with the perturber and wider resonances that they are 
captured into. 

We also calculate the orbital distributions of the parti- 
cles if their initial inclination angles to the perturber's 
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Fig. 12. — Orbital distribution of simulated test particles with initial inclination angles i G [20°, 30°] when the perturber with mass 
5 X 10^ Mq migrates to 0.3 pc. Other initial conditions are the same as those in Figure [D 



orbital plane i are in the range 20°-30°, 110°-120°, 
and 170°-180°, respectively. For the cases of initial 
i e [20°, 30°] and i € [110°, 120°], if the perturber mi- 
grates to a position Vp ~ 0.15 pc, the stars in the inner 
radii cannot be maintained on a disk (with 9f up to 50° 
or 120°, e.g., see Fig. [TT|) due to the nodal precession 
of the particles' orbital planes, which is consistent with 
our discussion in § 14.11 If the perturber has a smaller 
mass and migrates to a relatively outer position (e.g., 
Mp = 5 X 10^ Mq and = O.Spc), the angles 0f of 
the inner stars can be small (e.g., < 20°), but a warp- 
ness may be developed in the disk (the warpness is not 
indicated by observations so far). For the case of initial 
i e [20°, 30°] (see Fig.[T2l), both the eccentricities and in- 
clinations of the stars at the outer region can be excited 
up due to capture into mean-motion resonances, close en- 
counters with the perturber, and nodal precessions. But 
the eccentricities cannot be sufficiently excited up in the 
case of initial i € [110°, 120°] (e.g., with e/ < 0.6), be- 
cause our simulations show that they are less likely to 



be captured into mean-motion resonances compared to 
the case of i £ [20°, 30°] or that their mean-motion res- 
onances are more likely to be unstable and affected by 
close encounters with the perturber. A relatively higher 
perturber mass may increase the excited eccentricities 
due to close encounters, but it also induces larger nodal 
precession of stars at the inner radii. 

For the case of initial i G [170°, 180°], the inner stars 
may be maintained on a disk and on nearly circular or- 
bits. The angles 9f oi the stars with large initial semi- 
major axes (e.g., > 0.3 pc) can be excited up to 50°, 
and their eccentricities can also be excited up. But many 
excited particles have large a/ (e.g., > 0.5 pc) and the ex- 
cited eccentricity e/ of the particles at r ~ 0.15 — 0.5 pc 
arc not high enough (e.g., < 0.6) compared to the ob- 
servations shown in Figure 1131 In this case, both of the 
inclination and eccentricity excitations come mainly from 
close encounters with the perturber. 

If the perturber is a dark cluster (with a central 
IMBH), its orbital decay may be accompanied with mass 
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loss induced by the strong external tidal potential. Our 
simulations indicate that for the same initial perturber 
mass, a decrease in the perturber mass during its orbital 
decay may make the inner stellar disk more likely to be 
preserved. 

We also test the effects of different perturber- migration 
timescales. For longer migration timescales (e.g., 5Myr 
or lOMyr), the main results are generally not affected. 
Note that a very large lOMyr) migration timescale 
is not compatible with the estimated age of the young 
stars. For a substantially shorter migration timescale 
(e.g., r = IMyr), the inward- migration of the perturber 
is too fast for some particles to be captured into the 
mean-motion resonance and the excitation of the eccen- 
tricity and inclination becomes insufficient. 

4.3. Comparison with observations and discussions 

The numerical models in § 14.21 simulate the interac- 
tion between a disk of stars with nearly circular orbits 
and an inward-migrating perturber. Provided the per- 
turber's orbit is either nearly coplanar (in corotating di- 
rections) or essentially overhead, the young stars at the 
outer region generally attain high eccentricities (up to 
unity), while those at the inner region retain low eccen- 
tricities. In addition, the young stars at the outer region 
tend to have a large range of inclination angle relative to 
the disk, while those in the inner region roughly remain 
in the initial disk (see panels a and b in Figs. [7] and [8|). 

These features are consistent with the observed fea- 
tures of the orbits partly hsted in the items (a), (b) and 
(c) in § [T] In order to make a quantitative compari- 
son with observations, we plot the eccentricities of the 
young stars in the GC versus their observed distance to 
the central MBH in F igure[T3l (the data are adopted from 
iPaumard et al.l[2006h . This figure indeed indicates a cor- 
relation that the stars in the outlying regions of GC have 
high eccentricities (close to 1) while the inner ones only 
have moderate eccentricities. This apparent correlation 
is reproduced in our model as shown in panel (a) in Fig- 
ures [7] and [8] Note that the low eccentricity of the test 
particles in the inner region (with rf < 0.15 pc) preserved 
(panel a in Figs. [7] and [51), which appears inconsistent 
with the observational range of 0.2 — 0.4 (Fig. [T3)) . How- 
ever, the eccentricities of the young stars at the inner 
region may relax to the observed valu es due to the in- 
teract ions among the stars as shown in [Alexander et al.l 
()2007[ ). which are not considered in §|4l An N-body nu- 
merical simulation of the model proposed in this paper 
would provide a quantitative consistency check. Our cal- 
culations in § 14.21 show that the inclination of the stars 
in the outer region relative to the inner disk is less ex- 
cited for a nearly coplanar perturber orbit (only up to 
~ 50°) than for an overhead perturber orbit. The incli- 
nation angles of the outer stars are not well determined 
in observations, but if they have higher inclinations (e.g., 
~ 110° in item a in the §[T]), the perturber would be more 
likely to be on an overhead orbit. Precise determination 
of the inclination of the stars in the outer region will 
be important to distinguish the inclinations of the per- 
turber. 

In the in-spiraling-cluster scenario (briefly discussed in 
§[T]to account for the origin of the young stars), we note 
that, within a limited dispersion, the eccentricities of the 
tidally disrupted stars are usually comparable to that 
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Fig. 13. — The eccentricities of the young stars in the GC ver- 
sus their observed distance to the central MBH. The open cir- 
cles represent those stars having the measurements of the (three- 
dimensional) distance (r). The solid circles represent those stars 
only having the measurements of the projected distance (rpjoj). 
The inset shows several S-stars which have the eccentricity esti- 
mates. The figure shows a tendency that the outer stars have high 
eccentricities (close to 1) while the inner ones (with exclusion of 
the S-stars) only have moderate eccent ricities. The data used here 
are adopted from lPaumard et al.ll2006i . Using the Spearman corre- 
lation analysis, we find a quite strong correlation between the ec- 
centricity and the projected distance and the correlation coefficient 
between the eccentricity and the projected distance is = 0.58 
with Pflg = 8 X 10"'^. 

of the in-spiraling star cluster. The dominant compo- 
nent of the velocity vector of the dispersed stars is their 
co-moving motion with the cluster and the velocity of 
the stars relative to the center of the cluster is relatively 
small (Berukoff & Hansen 2006). Therefore, it is unlikely 
that the observed eccentricity distribution of the young 
stars shown in Figure [13] can be produced by a disrupted 
single in-spiraling young star cluster. In order to repro- 
duce the observed velocity distribution, the disruption 
of a second hypothetic al in-spiraling young star cluster 
may be required (e.g.. [Berukoff fc HansenI [2006f ) . Such 
a scenario seems contrived because the coeval nature of 
the young stars also requires these clusters to be formed 
simultaneously. 

The model we analyzed here is different from that two 
in-spiraling star cluster scenario. In our model, all the 
young stars are initially in a single disk and then the 
orbits of the young stars in the outer region of the disk 
are perturbed by an inward-migrating (dark) star cluster 
or an IMBH onto high eccentricity and high inclination 
orbits. 

Figures [71 [5] and [T^ (e) and (f) indicate one of the fea- 
tures of our model that because they are captured into 
mean-motion resonance, some test particles have nearly 
the same scmimajor axes with widely different eccentric- 
ities. The particles captured into the resonance are likely 
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to have close encounters with each other, which should 
modify the orbital distribution obtained here. Our pre- 
liminary N-body results show that some of the particles 
captured into the mean-motion resonance may then be 
scattered out of the resonance. Nevertheless, a substan- 
tial number of particles still attached to the resonance. 
The main conclusion derived in Figs. [7| and [5] are not 
expected to change significantly. The existence of res- 
onant stars, if can be extracted from the observational 
data, would not only provide an important check for our 
model but also strongly signify the possible presence of 
an IMBH in the GC. 

Our simulation results show that the inclination angles 
of the young stars at the outer region of the initial disk 
may be excited up to different values as the perturber 
migrating inward. The excited stars cannot form a co- 
herent secondary disk in this model. The formation of a 
secondary di sk is also challenging in th e in-spiraling clus- 
ter scenario (jBerukoff &: Hansenll2006( ) . In any event, ob- 
servational evidence of the existence of a secondary disk 
remains controversial. 

5. CONCLUSIONS 

We studied the dynamical evolution of a young stel- 
lar disk surrounding the MBH in the GC and perturbed 
by an inward migrating IMBH or a (dark) star cluster 
hosting an IMBH. Our numerical simulations show that 
the orbits of the young stars in the disk may be signifi- 
cantly modified from the outside in by the perturber. If 
the perturber has a low inclination angle to the initial 
disk, its migration from a region outside to ^ 0.15 pc 
may excite the eccentricities (up to unity) and the incli- 
nation angles of stars in the outer (0.2-0.5 pc) regions by 
capturing them into its mean-motion and secular reso- 
nances, forcing them to migrate with it, and/or closely 
encountered with them. Stars interior to this region pre- 
serve their initial coplanar structure. The overall dynam- 
ical distribution of the stars reproduces that observed on 
sub-parsec scale around the GC, i.e., an inner disk sur- 
rounded by a torus of highly eccentric and inclined stars. 



If the perturber migrates to ^ 0.3 pc on an orbit which 
is nearly perpendicular to the initial disk, the stars in 
the outer regions can still be excited to highly eccentric 
and inclined orbits through resonant capture and nodal 
precession induced by and close encounters with the per- 
turber. The inclinations can be excited to be in a larger 
range in the perpendicular case than those excited in the 
low-inclination case. In the perpendicular case, the stars 
in the inner region can also retain their low eccentricities 
and remain on the initial disk. These results reproduce 
many features of the observed orbital distribution of sub- 
parsec young stars. Note that the predicted position for 
the perturber here is an instantaneous value, and the 
perturber does not necessarily stop there but continues 
to migrate inward under the action of dynamical friction 
with time going on. 

Further measurements of the orbital parameters of the 
young stars in the GC (e.g., semimajor axes, eccentric- 
ities, and inclination angles) would provide important 
tests for the model proposed in this paper. It would 
also be useful to distinguish various scenarios proposed 
to account for the formation of the young stars. An im- 
portant confirmation for the model proposed in this pa- 
per is the discovery of a perturber [either an IMBH or a 
(dark) cluster hosting an IMBH] with a mass in the range 
3 X 10^ - 3 X 10'' Mq at r 0.15- 0.3 pc with an orbital 
plane which is either nearly parallel or perpendicular to 
the inner disk plane. 



We have benefited from discussions with Scott 
Tremaine. This work is supported by NASA (NAG5- 
12151, NNG06-GH45G), JPL (1270927), NSF(AST- 
0507424). Q.Y. acknowledges initial support from NASA 
through Hubble Fellowship grant HST-HF-01169.01-A 
awarded by the Space Telescope Science Institute, which 
is operated by the Association of Universities for Re- 
search in Astronomy. Inc.. for NASA, under contract 
NAS 5-26555. 



REFERENCES 



Alexander, R. D., Begolman, M. C, & Armitage, P. J. 2007, ApJ, 
654, 907 

Alexander, T. 1999, ApJ, 527, 835 
Alexander, T. 2005, Phys. Rep., 419, 65 

Beloborodov, A. M., Levin, Y., Eisenhauer, F., Gcnzel, R., 

Paumard, T., Gillessen, S., & Ott, T. 2006, ApJ, 648, 405 
Berukoff, S. J., & Hansen, B. M. S. 2006, ApJ, 650, 501 
Binney, J., & Tremaine, S. 1987, Galactic Dynamics (Princeton, 

NJ: Princeton University Press) 
Eisenhauer, F., et al. 2005, ApJ, 628, 246 
Genzel, R., et al. 2003, ApJ, 594, 812 
Gerhard, O. 2001, ApJ, 546, L39 
Chez, A. M., et al. 2003, ApJ, 586, L127 
Chez, A. M., et al. 2005, ApJ, 620, 744 
Goodman, J. 2003, MNRAS, 339, 957 
Giirkan, M. A., & Rasio, F. A. 2005, ApJ, 628, 236 
Hansen, B. M. S., & Milosavljevic, M. 2003, ApJ, 593, L77 
Hopman, C, & Alexander, T. 2006, ApJ, 645, 1152 
Kim, S. S., Figer, P. F., & Morrison, M. 2004, ApJ, 609, 123 
Levin, Y. 2007, MNRAS, 374, 515 
Levin, Y., & Beloborodov, A. 2003, ApJ, 590, L33 
Levin, Y., Wu, A., & Thommes, E. 2005, ApJ, 635, 341 



Nagasawa, M., Lin, D. N. C, & Ida, S. 2003, ApJ, 586, 1374 
Nayakshin, S. 2005, MNRAS, 359, 545 

Nayakshin, S., Dehnen, W., Cuadra, J., & Genzel, R. 2006, 

MNRAS, 366, 1410 
McMillan, S. L. W., & Portegies Zwart, S. F. 2003, ApJ, 596, 314 
Melia, F., & Falcke, H. 2001, ARA&A, 39, 309 
Misner, C. W., Thorne, K. S., & Wheeler, J. A. 1973, Gravitation, 

W. H. Freeman and Co., San Francisco 
Morris, M. 1993, ApJ, 408, 496 
Paumard, T., et al. 2006, ApJ, 643, 1011 

Perets, H. B., Hopman, C, & Alexander, T. 2007, ApJ, 656, 709 
Polnarev, S. & Rees, M. J. 1994, A&A, 283, 301 
Quinlan, G. D. 1996, NewA, 1, 35 

Ranch, K., & Tremaine, S. 1996, New Astronomy, 1, 149 

Sanders, R. H. 1992, Nature, 359, 131 

Schodel, R., et al. 2002, Nature, 419, 694 

Yu, Q., 2002, MNRAS, 331, 935 

Yu, Q., & Tremaine, S. 2001, AJ, 121, 1736 

Yu, Q., & Tremaine, S. 2003, ApJ, 599, 1129 



